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Abstract 

We consider the problem of modeling the dependence among many time series. We build high dimensional 
time- varying copula models by combining pair-copula constructions (PCC) with stochastic autorcgressive 
copula (SCAR) models to capture dependence that changes over time. We show how the estimation of this 
highly complex model can be broken down into the estimation of a sequence of bivariate SCAR models, 
which can be achieved by using the method of simulated maximum likelihood. Further, by restricting the 
conditional dependence parameter on higher cascades of the PCC to be constant, we can greatly reduce 
the number of parameters to be estimated without losing much flexibility. We study the performance 
of our estimation method by a large scale Monte Carlo simulation. An application to a large dataset of 
stock returns of all constituents of the Dax 30 illustrates the usefulness of the proposed model class. 

Keywords: Stock return dependence, time-varying copula, D-vincs, efficient importance sampling, 

sequential estimation 

JEL Classification: C15, C51, C58 

1. Introduction 



The modeling of multivariate distributions is an important task for risk management and asset al- 
location problems. Since modeling the conditional mean of financial assets is rather difficult, if not 
impossible, much research has foc used on modeling co nditional volatilities and dependen cies. The lit 



erature on multivariate GARCH ( Bauwens et al 



1994 . 



Yu and Mevei 



20061 ) and stochastic volatility models (jHarvev et al 



20061) offers many approaches to extend univariate volatility models to multivariate 
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settings. However, usually the resulting multivariate model makes the assumption of (conditional) multi- 
variate normality. Multivariate models based on copulas offer a popular alternative in that non-elliptical 
multivariate distributions can be constructed in a tractable and flexible way. The advantage of using cop- 
ulas to construct multivariate volatility models is that one is free with the choice of the marginal model, 
i.e. the univariate volatility model, and that it is possible to capture possibly asymmetric dependencies 
in the tails of different pairs of the distributions. In particular, lower tail dependen ce o ften needs to be 
accou nted for when measuring financial risks. Among many others. 



PattonI (|2009[ ) or 



Cherubini et al 



(j2004[ ). and references therein, give an overview of copula based models in financial applications. 

Two major drawbacks of most of the early applications of copula based models are that most studies 
focus on bivariate copulas only, limiting the potential for real world applications, and that the dependence 
parameter is assumed to be time-constant. This is in contrast to the empirically observed time-varying 
correlations. Each of these issues individually has been addressed in the literature in recent years. 
Larger dimensional copulas other than Gaussian or S tudent copulas h a ve be come available thro u gh th e 



introduction of hierarchi cal Archimedean copu las by 



factor copula models by 



Savu and Tredd |2010f) and 



Okhrin et al. 



Oh and PattonI (j201l[ ). or the class of pair copula constructions by 



2009) 



Aasetal 



(|2009l ). In particular the latter class, also called vine copula constructions, has become extremely popu- 



lar because of its flexibility and because of the possibility of estimating the large number of par a meter s 



sequentially. E x ample s of financial applications o f vine copula models are, e.g. 



ChoUete et al 



Difimann et al 



torn, or 



been introduced by 



Patton 



(120091) 



Brechmann and Czadd (|201l[ ). Copulas with time- varying parameters have 



20061 ) to model changing exchange rate dependencies. Since then a number 



of stud ies have proposed different ways to specify time- v arying copulas. F o r exam ple 



(j2004l ) test for structural breaks in copula parameters 



point te sts to identify intervals of constant dependence, 



Giacomini et al 



Garcia and Tsafack 



3l 



Dias and Embrechts 



( 20091) u se a sequence of break- 



2011) and 



Stober and Czado 



(|2011al ) use a regime-switching model for changing dependencies, iHafner and Reznikoval (|2010l ) treat the 



copula parameter as a sm ooth function of time and estim ate it by local maximum likelihood, whereas 



Hafner and Manneij (|201l[ ) and I Almeida and Czadd (|201lf ) proposed a model where the copula parameter 



is a transformation of a latent Gaussian aut oregressive process of order on e. An overview and compari- 



son of time- varying copula models is given in 
for time- va rying parameters in larger dim ensions. 



Manner and Reznikova J 20111 ) . Only very few papers allow 



Chollete et al 



copula and 



of the DCC model by 



Engle 



20091 ) estimate a regime-switching vine 



Heinen and Valdesogo (j2008a[) allow the parameter of a vine copula to be driven by a variation 



20021) 
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The contribution of th is pa per is to extend the s tocha stic autoregressive copula (SCAR) model by 



Hafner and Manner! ( 2011 ) and 



Almeida and Czadd (|201ll) to practically relevant dimensions using D- 



vines. We discuss how the proposed model can be estimated sequentially using simulated maximum 
likelihood estimation. We also address how a number of conditional copulas can be restricted to be time- 
constant or independence copulas without restricting the flexibility of the model too much. Furthermore, 
we perform a large scale Monte Carlo study investigating the behavior of the proposed estimator and 
find that it shows an acceptable performance. In our empirical study we apply our model to the returns 
of 29 constituents of the German DAX 30 and find that the model performs quite well. 

The remainder of the paper is structured as follows. The next section introduces copulas in general, 
SCAR models, D-vines copulas and shows how they can be combined to obtain the ficxible class of D-vine 
SCAR models. Section [3] treats the estimation of the proposed model. Section 2] presents the results of 
our simulation study. Section [5] contains the empirical application and Section [6] gives conclusions and 
outlines further research. 

2. D-vine based SCAR models 

We are interested in modeling the joint (conditional) distribution of a d-dimensional time series 
Vt = (yi.t, Vd.t) for t = 1, T. We assume that each variable yi^t for i = 1, d follows an ARMA(p,q)- 
GARCH(1,1) process, i.e. 

p g 

yi,t = ft,0 + X] l^hjyi,t~3 + X! ^i,k'^ht-k£i,t-k + <^t,t£i,t (1) 
3 = 1 k=l 

with 

The usual stationarity conditions arc assumed to hold. Denote the joint distribution of the standardized 
innovations Si^t by G(ei_t, £d,t) and let their marginal distributions be Fi(ei^t), Fd{sd,t), respectively. 
Then by Sklar's theorem there exists a copula G such that 

C?(ei,t, ed,t) = C{F,{ei.t). Fd{£d,t))- (2) 

Since all the marginal behavior is captured by the marginal distribution, the copula captures the com- 
plete contemporaneous dependence of the distribution. Let Ui_t ~ Fi{ei^t) be the innovations transformed 
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to U(0, 1) random variable and define Ut := {ui^t, Ud.t)- 

In the remainder of this section we first show how time-varying dependence can be incorporated in 
bivariate copula models. Then we discuss how these models can be extended to arbitrary dimensions, for 
which we need the notion of vines. 

2.1. Bivariate SCAR Copula Models 

For now, consider the bivariate time series process {ui f,Uj t) for t = 1,...,T. We assume that its 
distribution is given by 

(uu,w,.t)~C(-,-;C) (3) 

with 91^ e e the time -varying parameter of the copula C. In order to be able to compare copula parameters 
that have different domains, the copula can equivalently be parameterized in terms of Kendall's r e 
(—1,1). This follows from the fact that for all copulas we consider there exists a one-to-one relationship 
between the copula parameter and Kendall's r, which we express by 91-' ~ r{Tl-'). We assume that Tj-* is 
driven by the latent Gausian AR(1) process given by 



(4) 



where z^-.t are independent standard normal innovations. We further assume \4>ij\ < 1 for stationarity 
and (T.jj > for identification. Due to the fact that Aj"* takes values on the real line we apply the inverse 
Fisher transform to map it into (—1, 1), the domain of Tj-': 



exp{2\. 



(5) 



2.2. D-vine Distributions and Copulas 



only the class of bivariate copulas (e.g. 



Joe 


1997 


and 


Nelsen 


2006) 



20061 ) was flexible enough to accommodate 



asymmetric and/or tail dependence without placing unrealistic restrictions on the dependence structure. 
Recently pair copula constructions (FCC) are found to be very useful to construct flexible multivariate 
copulas. Here a multivariate copula is built up with bivariate copula terms modeling unconditional and 
conditional dependencies. The first such constru ction was proposed i n 



sign ificantly extended to more ge neral settings in 



and 



Bedford and Cookd (|2002|) 



Joe 



(11996^. It was subsequently 



Bedford and Cookd mm 



Kurowicka and Cookd (|2006l ). They called the resulting distributions regular (R) vines and explored 
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them for the ease of Gaussian pair copulas. The backbone is a graphical representation in form of a 
sequence of linked trees identifying the indices which make up the multivariate copula. In particular, they 
proved that the specification of the corresponding pair copula densities make up a valid multivariate copula 
density. Further properti es, estimation, model sele ction methods and their use in complex modeling 



situations can be found in 



Aas at al. 



Kurowicka and Jod ( 2011 ) 



(120091 ) recognized the potential of this construction for statistical inference and developed 



a sequential estim ation (SE) procedur e , whi ch can be used as starting values for maximum likelihood 



estimation (MLE). [Bedford and Cookd (j2002l) identified two interesting subclasses of regular vines called 
D-vines and canonical (C)-vines. In the case of D-vines the sequence of vine trees consist of pair trees, 
while for C-vines they are starlike with a central node. This shows that C-vines arc more useful for data 
situations where the importance of the variables can be ordered. This is not the case for the application we 
will present later; therefore we concentrate on D-vines. However, we would like to note that multivariate 
SCAR models can also be constructed based on C-vines and more generally on R- vines. 

Notably, C- and D-vines can be introduced from first principles (e.g. 
(Xi, Xd) be a set of variables with joint distribution F and density /, respectively. Consider the 
recursive decomposition 



2010f) . For this let 



/(Xi, ...,Xd) 



Y[ fi^klxi, . . .,Xk-l) X f{xi). 



(6) 



k=2 



Here F{-\-) and later /(•!•) denote conditional cdf's and densities, respectively. As a second ingredient 
we utilize Sklar's theorem for dimension d — 2 to express the conditional density of Xi given X2 = X2 as 



/(a;i|a;2) = ci2(Fi(a;i), i^2(a;2)) x fi{xi), 



(7) 



where C12 denotes an arbitrary bivariate copula density. For distinct indices i, j, ii, ■ • ■ , ik with i < j and 
ii < • • • < ifc we now introduce the abbreviation 



\D ■= C,^j\DiFix^\TlD), F{xj\xD)), 



(8) 



where D := {zi, • • ■ , ik} and x/j :— (x^^ , 



J. Using ([7]) for the conditional distribution of {Xi,Xk) 
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given X2 = X2, ■ ■ ■ Xk-i = Xk-i we can express f{xk\xi, • • • , Xk-i) recursively as 



f{xk\xi,...,Xk-l) = Ci^k\2:k-lX f{xk\x2,---,Xk-l) 

k-2 



[H Cs,k\s+l:k-l] X C(k-l),k X fk{xk), 



(9) 



s=l 



where r : s := (r, r + 1, . . . , s) for integers r and s with r < s. Using ([9]) in ([6|) and s = i,k = i + j it 
foUows that 



d—l d—j d 

f{xi,...,Xd) = [[[Ylc,^^+]\^+l■.^+J-l]■[Y[fkixk)] 

j = l 4=1 k=l 



(10) 



If the marginal distribution of Xk are uniform for all /c = 1, • • ■ ,d, then we call the corresponding density 
in (jlOp a D-vine copula density and the corresponding distribution function a D-vine copula. 
For illustration we consider a five dimensional D-vine, its density then given by 



/(Xi,---,X5) = fk{xk)] ■ C12 ■ C23 ■ C34 



k=l 



X C45 • Ci3|2 • C24|3 ■ C35|4 • Ci4|23 • C25|34 • Cl5|234> 



(11) 



with the corresponding vine tree representation identifying the utilized indices given in Figure [T] In 
particular the indices in Tree Ti indicate the unconditional pair copulas, while Trees T2 to correspond 
to conditional pair copulas, where the set of conditioning variables has size 1 to 3, respectively. 
If Ci^i^j\i^i,i^j_i models the dependence between the rv's F{Xi\Ki+i:i+j-i) and 

|xi_|_i:i_|_j_i ) we implicitly assume that the copula density (•, •) does not depend on 

the conditioning variables jCj+i-j+j^ i other than through the arguments |xi-)_i:i_|_j_i) and F{Xi^j |xi-)_i:i_|_j_i). 



This is a common assumption and 



Haffetal 



(|201C1I) call this a simplified vine. They showed that this 
restriction is not severe by examining several examples. 

In the D-vine representation given in (1101) we a lso need a fast recursive way to compute conditional 



cdf's which enter as arguments. For this 



Jod (119961 ) showed that for w G and := D\v 



(12) 
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12 



23 



34 



45 



5^ Ti 



12 ; 1 23 j 1 34 J ^ 45 , 

1312 2413 3514 




,1312 



14123 
14123^ 




2413, 




151234 



25134 
25^ 




3514, 




Figure 1: A D-vinc tree representation for d = 5. 



For the special case of D = {t;} it follows that 



In the case of uniform margins uj — Fj{xj), for a parametric copula cdf Cjy(uj,u,j) = Cjy{uj,Uv]Ojy) 
this further simplifies to 



h{uj\uy,6j^) : 
With this notation we can express F(xj\xo) as 



F{xj\xD) = h{F{xj\xD^J\F{xy\xD^J,6jv\D.J- 



This allows the recursive determination of the likelihood corresponding to (jlOp . Furt hermore, the inv erse 



of the fe-functions is used to facilitate sampling from D- and C-vines (see for example 



Kurowicka and Cooke 


2007) 


Stober and CzadJ 


2011h 


). 



Aas et al 



2009 and 



20071) . They are also used for sampling from the more general R-vine model (see 
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2.3. D-vine based multivariate SCAR models 

We now combine bivariate SCAR models and D-vines to formulate a multivariate D-vine SCAR model. 
For this we use a bivariate SCAR copula model as the pair copula model in a D-vinc copula. This gives 
rise to the following definition of a D-vine SCAR copula density 

d-ld-j 

j=i i=\ 

where := + l : i+j — 1 and Qt '■= {^t'*'"'^ j = l,---,<i — l,i = 1, ■ ■ ■ , d ~ j} is the 

time-varying copula parameter vector. Here c(-, •; is the bivariate copula density corresponding to 

the bivariate SCAR copula given in ([3]), where satisfies 

e^-''^=r-(r,'('^^') = r(V'(Al(^'^))) (15) 

for the latent Gaussian AR(1) process A^'*'"''' with 

Here z'^^'''^ are independent standard normal innovations for j = 1, • ■ • , d; i = 1, • • ■ , d — j. As for the 
bivariate case, we assume < 1 and cr^i j) > for stationarity and identification, respectively. The 

bivariate copula family corresponding to l{i,j) can be chosen arbitrarily and independent of any other 
index lir, s). 

The copula in p4)) can be used in ([2]) to specify the joint distribution of the innovations in ([1]). 

3. Parameter estimation in D-vine SCAR models 

We are interested in estimating the parameters of both the marginal models and the stochastic copula 
models. The joint density of our model is given by the product of the marginal and the copula densities 

5(ei,t,...,£d,t) = c(Fj(ei,t), ...,Fd{ed,t)) ■ Mei.t) ■ •■• • fd{£cLt), 

where g, c and / denote the densities of the joint distribution, the copula and the marginal distributions, 
respectively. Taking logarithms, we can see that the joint log-likelihood is the sum the marginal and 
the copula log-likelihood function. For estimation we ultilize a two-step approach common in copula 
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based models. In this approach first the marginal parameters are estimated 
residuals are formed. These are transfo rmed using either a parametric (see 



separately and standardized 



Joe 



probability integral transformation (see 



Genest et al. 



20051 ) or nonparametric 



19951 ) to get an independent sample from a multi- 



variate copula. These transformations do not change the dependence structure among the standardized 
residuals. This approach allows us to perform the estimation of the marginal and copula parameters 
separately. If the marginal models are chosen carefully, as we will do. then a parametric probability 
transformation is a good approximation to the t rue copula data n^.f = Fi[ei^t)- Problems only occur if 



the marginal models are grossly misspecified (see 



Kim et al 



20071) 



Furthermore, we saw above that the density of a D-vine copula is the product of bivariatc (conditional) 
copulas. Therefore, instead of estimating all copula parameters of our model in one step, which is 
computationally infeasible due to the large number of parameters, we are able to estimate the copula 
parameters sequentially. 

In Section 13.11 the estimation of bivariate SCAR copula models by simulated maximum likelihood 
(SML) using efficient importance sampling (EIS) is reviewed. Section [3.21 presents the sequential estima- 
tion of vine copula models and in Section 13.31 we discuss how the sequential estimation of D-vine SCAR 
copula models can be achieved. 

3.1. Estimation of bivariate SCAR copula models 

For the moment, we are interested in estimating the copula parameter vector lj := (/i, ct). For 
notational convenience we decided to drop the indices i and j whenever no ambiguity arises. Denote 
Ui ~ {ui^t}t=i, Uj = {uj^t^J^i and A = {\t}J=i and let f{ui,Uj,A;uj) be the joint density of the 
observable variables {ui,Uj) and the latent process A. Then the likelihood function of the parameter 
vector cl) can be obtained by integrating the latent process A out of the joint likelihood, 



(17) 



We can alternatively write this as a product of conditional densities 



■J + — 1 



(18) 



This is a T-dimcnsional integral that cannot be solved by analytical or numerical means. It can. 
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however, be solved efBci ently by Monte Car l o inte gration using a technique called efficient importance 



sampling introduced by iRichard and Zhang (l2007|) . The idea is to make use of an auxiliary sampler 
m{Xt; At_i, at) that utilizes the information on the latent process contained in the observable data. Note 
that it depends on the auxiliary parameter vector at = (ai,t, a2.t)- Multiplying and dividing by m(-), the 



likelihood can then be rewritten as 



^ '"/(ui,t,-Uj,t, At|At_i,a;)"' ^ 



W_■m{\t]\t-l,at)d^.. 



m{\t]\t-i,at) 

Drawing N trajectories A^'-' from the importance sampleiQ the likelihood can be estimated by 



(19) 



N / T 



i=l \t=l 



m(A«;A«i,a,) 



(20) 



This leaves the exact choice of the importance sampler m(-) to be determined, which ideally should 
provide a good match between the numerator and the denominator of (j20p in order to minimize the 
variance of the likelihood function. It is chosen to be 



m{Xt;Xt-i,at) 



fc(Af, At-i; ot) 
x(At_i;at) 



(21) 



where 



is the normalizing constant of the auxiliary density kernel fc(-). Furthermore, the choice 



x(At_i;at)= / k{Xt,\t-i;at)d\t 



k{Xt,Xt-i;at) = p(At|At_i, a;)C(At, at), 



with p(Af I At_i, w) the conditional density of At given Xt-i and C,{Xt,at) = exp(ai.tAt +a2.tXt) turns out 
to simplify the problem considerably. Noting that /(w^^i, Wj^i, At|At_i, u;) — c{ui^f,Uj^tT ^t)p{^t\^t-i,'^), 
the likelihood expression (|19p can be rewritten as 



L{uj;u,,uj) = / n 



T r 



c(uj,t,Wi,t; A()x(At;a4+i) 



exp(ai,tAt + a2.tXt) 



Y\_m{Xt; Xt-i, at)dA, 



(22) 



good choice for N is about 100. 
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where we have used the fact that x(-) can be transferred back one period, because it does not depend 
on At. Defining xi^Ti Qt+i) = 1 and given a set of trajectories A^*) for i — 1, . . . , N, minimizing the 
samphng variance of the quotient in the hkehhood function is equivalent to solving the following linear 
least squares problem for each period t = T, . . . , 1, 

log c(ui,t, Uj,t; A^^) +logx(At;at+i) = ct +ai,fA|*'' + ai,t[A^*'']^ + Vt'^ ■ (23) 

This problem can be solved by OLS with ct the regression intercept and yy^ the error term. Then the 
procedure works as follows. First, draw N trajectories A^*) from p(At|Af_i, w) and estimate the auxiliary 
parameters at for t = T, . . . , 1 by solving (|23p . Next, draw N trajectories A^'^ from the importance 
sampler m(At; Af-i, at) and re-estimate the auxiliary parameters {dj}^]^. Iterate this procedure until 
convergence of {dt}J^^ and use N draws from the importance sampler to estimate the likelihood function 
(PU)) . This likelihood function can then be maximized to obtain parameter estimates a>. Note that 
throughout the same random numbers have to be used in order to ensure convergence of {at}J^i and 
smoothness of the likelihood function. 

Although the parameter vector a; driving the latent process is of some interest, ultimately one wishes 
to get estimates of the latent process A and transformations thereof. In particular, we are interested 
in estimating Tt = 'ip{Xt) for t = 1, . . . ,T, where ip{-) denotes the inverse Fisher transform given in ([5]). 
Smoothed estimates of ip{Xt) given the entire history of the observable information Ui and Uj can be 
computed as 



E[4}{\t)\Ui,Uj\ ^ ^ — ^ _ ^^^^ . (24) 



/- 0(At)/(M,,Mj-, A;a;)rf A 

Note that the denominator in corresponds to the likelihood function and both integrals can be 
estimated using draws from the importance sampler m(At; At_i, at). Filtered esti mates of il'iM) given in- 



forma tion until time t—1 can be computed in a similar way and details are given in 



Liesenfeld and Richard 



mm 



3.2. Sequential estimation of D-vine copula parameters 

The form of the D-vine density given in (|10p allows for a sequential p arameter estimati on approach 
starting from the first tree until t he last tree. This w as first proposed by 
and shown in detail for C-vines in 



Aas et al. 



Czadoetal 



( 2009 ) for D-vines 



()2011l ). First estimate the parameters corresponding to 



the pair-copulas in the first tree using any method you prefer. For the copula parameters identified in the 
second tree, one first has to transform the data with the h function in (fT3)) required for the appropriate 
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conditional cdf using estimated parameters to determine pseudo realizations needed in the second tree. 
Using these pseudo observations the parameters in the second tree are estimated, the pseudo data is 
again transformed using the h function and so on. 

For example we want to estimate the parameters of copula Ci3|2. First transform the observations 
{ui,t,U2,t,U3^t,t = I,-- - ,n} to Ui\2,t ■= ^(ui,t|"2,t,^i2) and U3|2,t := /i(u3,t|w2,t, ^23), where §12 and §23 
are the estimated parameters in the first tree. Now estimate 0i3|2 based on {wi|2,t, U3|2,t; t — 1, • • • , n}. 
Continue sequentially until all copula parameters of all trees are estimated. For trees Ti with i > 2 
rec ursiv e appli cations of the h functions is needed. Asymptotic normality of the SE has been established 



by 



Hafn (|2010f ). However, the asymptotic covariance of the parameter estimates is very complex and 



one has to resort to bootstrapping to es ti mate the standard err o rs. SE is often used in large di men- 



sional problems, e.g. 



Mendes et al. 



Brechmann and Czado 



torn . 



Heinen and Valdesogol (|2008bf ) 



Brechmann et al 



( 20111) and 



201 ih . A joint MLE of all parameters in (fTpl) requires 



mizat ion. Therefore SE's are often used as starting values as, e.g., in 



(I2OIII) 



Aas et al 



high dime nsional opti- 



(I2OO9I) and 



Czado et al 



A final issue to be discussed is model selection for D-vine copula models. For non Gaussian pair 
copulas, permutations of the ordering of the variables give different D-vine copulas. In fact, there are 
d\/2 different D-vine copulas when a common bivariate copula is used as pair copula type. Often the 
bivariate Clayton, Gumbel, Gauss, t, Joe and Frank copula families are utilized as choices for pair copula 
terms. However, in this study we restrict the attention to the Gauss, Gumbel, Clayton and rotated 
versions thereof. 



3.3. Estimation of D-vine SCAR models 

In principle, estimation of the D-vine SCAR model given in (|14p works the same way as for static 
D-vinc copulas. There are, however, two important differences. First of all, given that the bivariate 
SCAR models in the first tree have been estimated, it is not possible to apply the h function given in 
(jl3p directly to obtain the pseudo observations that are needed to obtain the parameters on the second 
tree. The reason is that one only obtains parameter estimates of the hyper-parameters (p, (p, cr), but not 



of the latent (time-varying) copula parameters Ot- We do, however, have N simulated trajectories 9 



^From a practical perspective, the recent R package CDVine of lSchepsmeier and BrechmannI 1I2OIII ) provides easy to use 
random number generation, and both SE and MLE fitting algorithms for C- and D-vines. 
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from the importance sampler. With these we can calculate the pseudo observations by 

(25) 



N 

i 



N 

i=l 



where we suppress the dependence of ^ on the variable indices j and v for notational reasonso The second 
difference to the static D-vine copula model is that one-step estimation by MLE is computationally not 
feasible for the time- varying model, because each bivariate likelihood function needs to be computed by 
simulation. 

For a d dimensional dataset the D-vine SCAR copula one has to estimate 3d{d + l)/2 parameters. 
Fortunately, we can reduce the number of para meters to b e estim ated by placing a number of restrictions. 
Similar to tail properties of D-vines studied in 



Joeetal 



(|2010l ) the choice of time- varying pair copulas 



in the first tree propagates to the whole distribution, in particular all pairs of variables have an induced 
time-varying Kendall's tau. We expect estimation errors to increase for parameters as the corresponding 
tree increases because of the sequential nature of the estimation procedure. Therefore we allow for D-vine 
SCAR copula models where the pair copulas are time- varying only in lower trees, while the pair copulas 
are time-constant for higher trees. Such models will also be investigated in our simulation study. 

A second useful restriction is to allow for the possibility of truncating the D-vine copula, which means 
that we set all pair copulas beyond a certain tree equal to the independence copula. This is empirically 
justified, since the dependence in the lower trees seems to capture most of the overall dependence in 
the data and the conditional dependence in higher trees is hardly visible. Note that this also allows 
the estimation of our model in arbitrarily large dimensions, as we will only need to estimate (bivariate) 
models up t o a certain dimens i on an d can truncate the model thereafter. For static models this has been 
followed by 



Brechmann et al 



(|201l[ ) and includes tests at which level to truncate. 
In order to decide which copula family to use and whether to use time-varying, time-constant or 
independence copulas at certain levels we compare the Bayesian Information Criterion (BIG) for all 
competing models. We decided for this information criterion, because it favors parsimonious models. 
Given the high fiexibility of the D-vinc SCAR model and the difficulty to estimate the parameters at 
higher level, we believe that parsimony is crucial. 



^Alternatively, we could calculate the pseudo realizations using the smoothed estimates of the latent dependence pa- 
rameter using II 241 1. However, averaging over the nonlinear transformation h seems more reasonable than applying the 
transformation to the (weighted) average. 
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3.4- Computational issues 

Estimation of a bivariate SCAR model by simulated maximum likelihood programmed in MATLAB 
can take up to several minutes on a normal computer. As we consider a Monte Carlo simulation and 
an application with 29 variables this is much too slow. However, the problem at hand offers itself to 
parallel computing. On each tree one has to estimate a large number of bivariate models independently 
of each other. Therefore, given a sufficient number of processing cores we can estimate all models on one 
tree at the same time and proceed to the next tree once all models are estimated. Depending on the 
dimension of the problem, this can lead to immense increases in computing speed. The computations for 
the Monte Carlo study in the next section and our empirical application were performed on a Linux cluster 
computer computer with 32 processing cores (Quad Core AMD Opteron, 2.6Mhz), The most demanding 
computational task, the estimation of the log-likelihood function by EIS, is implemented in CH — h, which 
resulted in our code being about 20-30 times faster compared to MATLAB code. The maximization of 
the likelihood and the parallel computation within levels is implemented in R (version 2.12.1) by using 
the optim function and the multicore library. 

4. Simulation study 

To investigate the performance of the sequential estimation utilizing EIS for D-vine SCAR copula 
models we conducted an extensive simulation studyQ For this we chose a four dimensional setup. We 
simulated 1000 four dimensional SCAR D-vine copula data sets of length 1000 under several scenarios. 
Then we estimated the relative bias and relative MSE for the parameters of the latent AR(1) processes 
corresponding to each of the six bivariate copula terms of the D-vine copula. Standard errors are also 
estimated and given in parentheses in the tables. 

We expect that the stationary latent AR(1) signal-to- noise ratio given by sn :— ^f^^_^2-^-i/2 a-nd the 
stationary variance avar := i_^2 will influence the performance. Therefore we include these quantities 
as well. The copula families we consider arc the Gaussian, Clayton and Gumbcl copulas. 

In the following we present the small sample performance results for scenarios involving a single 
bivariate copula family for all pair copula terms and a common time-varying parameter structure for all 
pair copula terms (Section 14. ip . a single bivariate copula family for all pair copula terms and a common 
time- varying parameter structure for only pair copula terms in the first tree (Section 14. 2p . and different 

*For brevity we only present a part of the results and the outcomes for further parameter and model constellations are 
available from the authors upon request. 



14 



time- varying structures for all or only the first tree pair copulas and/or different copula families (Section 
14. 3[) . In Scction r4.4l we draw conclusions from our simulation study relevant for the subsequent application 
of our model. 

4-.1. Common copula families and common time-varying structure for all pair copulas 

The results in Table [1] show satisfactory results for all latent AR(1) parameters of unconditional pair 
copula terms, i.e. the terms with indices 12,23,34. The mean parameters fi and the standard error 
parameters a of the pair copulas in trees 2 to 3 are also well estimated in scenarios where the asymptotic 
AR(1) signal-to-noise ratio sn is large. The estimation of the persistence parameter is less effected by 
the value of the signal-to-noise ratio. By the sequential nature of the estimation procedure we expect the 
performance to deteriorate for parameters with increasing conditioning set. This behavior is also visible. 
In particular, the estimation on the second tree is still reasonably precise for most cases, but on the third 
level the results become significantly worse. Comparing the results across different copula families one 
can generally state that for the Gaussian copula the results are best, closely followed by the Gumbel 
copula. For the Clayton copula the estimation is most imprecise both in terms of relative bias and MSE. 

4-.2. Common copula families and common time-varying structure for pair copula terms in the first tree 
only 

The small sample performance for scenarios where only the copula parameter for pair copulas in the 
first tree is time- varying, reported in Table [51 is quite satisfactory for all models except for the mean 
parameters fi in the third tree. This was to be expected by the results from the previous section. Note, 
however, that mean dependence, measured by the parameter fj., in the second and third tree is lower 
than in the above scenario. We can therefore conclude that the precision in estimating /i decreases when 
the degree of dependence decreases. When the overall dependence is comparable, as is the case for the 
pair 13|2 in the second tree, we see that the parameter fi is estimated much more precisely than for the 
time-varying case. Furthermore, note that the average relative bias on the higher trees is positive for the 
copulas with lower fi parameter. 

4-3. Different time-varying structure and mixed or common pair copula families 

Since many scenarios are possible in this setup, we restrict ourselves to a three dimensional setup 
assuming all pair copulas time-varying with common pair copula family or different pair copula families. 
In particular, we vary the persistence parameter, which is now lower for one copula on the first tree and 
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Index 


True values 


Average relative Bias 




Average relative MSE 


sn 


avar 




^ (j) (7 


^ 4> 


a 


/i. (T 







Copula families: Gaussian Gaussian Gaussian Gaussian Gaussian Gaussian 



12 


.50 


95 


.15 


.0298(.0071) - 


0116(.0008) - 


0761(.0033) 


.0503( 


.0041) 


.0008( 


.0001) 


.0166( 


.0009) 


1.04 


0.23 


23 


.50 


95 


.15 


.0361(.0075) - 


0123(.0009) - 


0743(.0035) 


.0565( 


.0048) 


.0009( 


.0001) 


.0179( 


.0011) 


1.04 


0.23 


34 


.50 


95 


.15 


.0335(.0077) - 


0108(.0008) - 


0800(.0034) 


.0603( 


.0059) 


.0008( 


.0001) 


.0176( 


.0008) 


1.04 


0.23 


13|2 


.50 


95 


.15 


-.1036(.0051) - 


0061(.0006) - 


2112(.0036) 


.0362( 


.0016) 


.0004( 


.0000) 


.0573( 


.0016) 


1.04 


0.23 


24 3 


.50 


95 


.15 


-.0997(.0054) - 


0064(.0006) - 


2154(.0034) 


.0385( 


.0016) 


.0003( 


.0000) 


.0581( 


.0016) 


1.04 


0.23 


14|23 


.50 


95 


.15 


-.4090(.0042) - 


0215(.0011) - 


3230(.0050) 


.1846( 


.0035) 


.0017( 


.0001) 


.1286( 


.0031) 


1.04 


0.23 


Copula families 


: Clayton 


Clayton Clayton 


Clayton Clayt 


on Clayton 


















12 


.50 


95 


.15 


.0580(.0084) - 


0092(.0011) - 


0372(.0040) 


.0682( 


.0103) 


.0012( 


.0001) 


.0162( 


.0011) 


1.04 


0.23 


23 


.50 


95 


.15 


.0743(.0085) - 


0096(.0011) - 


0384(.0041) 


.0723( 


.0095) 


.0013( 


.0002) 


.0172( 


.0014) 


1.04 


0.23 


34 


.50 


95 


.15 


.0592(.0084) - 


0086(.0010) - 


0394(.0040) 


.0684( 


.0094) 


.0011( 


.0001) 


.0165( 


.0012) 


1.04 


0.23 


13|2 


.50 


95 


.15 


-.1915(.0047) - 


0110(.0008) - 


2585(.0040) 


.0572( 


.0022) 


.0007( 


.0001) 


.0818( 


.0022) 


1.04 


0.23 


24 3 


.50 


95 


.15 


-.1956(.0048) - 


0115(.0007) - 


2615(.0038) 


.0596( 


.0022) 


.0006( 


.0000) 


.0816( 


.0021) 


1.04 


0.23 


14|23 


.50 


95 


.15 


-.5676(.0036) - 


0379(.0019) - 


4283(.0059) 


.3338( 


.0039) 


.0048( 


.0007) 


.2148( 


.0046) 


1.04 


0.23 


Copula families 


: Gumbcl Gumbcl Gumbcl Gumbcl Gumbcl Gumbcl 


















12 


.50 


95 


.15 


.0210(.0063) - 


0036(.0006) - 


0220(.0036) 


.0393( 


.0024) 


.0004( 


.0001) 


.0133( 


.0007) 


1.04 


0.23 


23 


.50 


95 


.15 


.0258(.0064) - 


0045(.0006) - 


0238(.0039) 


.0417( 


.0025) 


.0004( 


.0001) 


.0154( 


.0015) 


1.04 


0.23 


34 


.50 


95 


.15 


.0142(.0064) - 


0038(.0006) - 


0267(.0037) 


.0409( 


.0032) 


.0004( 


.0000) 


.0143( 


.0006) 


1.04 


0.23 


13|2 


.50 


95 


.15 


-.1196(.0053) - 


0045(.0006) - 


2187(.0040) 


.0424( 


.0024) 


.0004( 


.0000) 


.0635( 


.0018) 


1.04 


0.23 


24 3 


.50 


95 


.15 


-.1168(.0052) - 


0050(.0006) - 


2217(.0038) 


.0405( 


.0018) 


.0004( 


.0000) 


.0633( 


.0018) 


1.04 


0.23 


14|23 


.50 


95 


.15 


-.4553(.0040) - 


0220(.0012) - 


3703(.0052) 


.2234( 


.0036) 


.0020( 


.0002) 


.1643( 


.0037) 


1.04 


0.23 


Copula families 


: Gaussian Gaussian Gaussian Gaussian Gaussian Gaussian 
















12 


.50 


95 


.05 


-.0023(.0022) - 


0177(.0017) 


0410(.0112) 


.0050( 


.0003) 


.0033( 


.0005) 


.1325( 


.0097) 


3.12 


0.03 


23 


.50 


95 


.05 


.0014(.0023) - 


0169(.0020) 


0459(.0107) 


.0055( 


.0003) 


.0043( 


.0021) 


.1227( 


.0081) 


3.12 


0.03 


34 


.50 


95 


.05 


.0029(.0023) - 


0164(.0016) 


0506(.0109) 


.0057( 


.0003) 


.0030( 


.0005) 


.1274( 


.0105) 


3.12 


0.03 


13|2 


.50 


95 


.05 


.0057(.0023) - 


0192(.0022) 


0704(.0111) 


.0058( 


.0003) 


.0053( 


.0031) 


.1354( 


.0081) 


3.12 


0.03 


24 3 


.50 


95 


.05 


.0013(.0022) - 


0214(.0020) 


0692(.0123) 


.0053( 


.0003) 


.0047( 


.0010) 


.1636( 


.0144) 


3.12 


0.03 


14|23 


.50 


95 


.05 


-.0740(.0023) - 


0310(.0029) 


0417(.0138) 


.0109( 


.0006) 


.0100( 


.0036) 


.2024( 


.0154) 


3.12 


0.03 


Copula families 


: Clayton 


Clayton Clayton 


Clayton Clayt 


on Clayton 


















12 


.50 


95 


.05 


-.0028(.0022) - 


0163(.0016) 


0465(.0105) 


.0049{ 


.0003) 


.0028( 


.0005) 


.1161( 


.0086) 


3.12 


0.03 


23 


.50 


95 


.05 


.0011(.0022) - 


0155(.0014) 


0586(.0104) 


.0052( 


.0002) 


.0021( 


.0003) 


.1170( 


.0091) 


3.12 


0.03 


34 


.50 


95 


.05 


.0030(.0023) - 


0141(.0012) 


0531(.0098) 


.0054( 


.0003) 


.0017( 


.0002) 


.1029( 


.0067) 


3.12 


0.03 


13|2 


.50 


95 


.05 


-.0365(.0021) - 


0162(.0013) - 


0321(.0099) 


.0060( 


.0003) 


.0021( 


.0002) 


.1037( 


.0061) 


3.12 


0.03 


24 3 


.50 


95 


.05 


-.0415(.0021) - 


0157(.0014) - 


0528(.0103) 


.0061( 


.0003) 


.0024( 


.0004) 


.1134( 


.0069) 


3.12 


0.03 


14|23 


.50 


95 


.05 


-.2252(.0019) - 


0727(.0043) 


1145(.0175) 


.0544( 


.0009) 


.0249( 


.0043) 


.3317( 


.0199) 


3.12 


0.03 


Copula families 


: Gumbcl Gumbcl Gumbcl Gumbcl Gumbcl Gumbcl 


















12 


.50 


95 


.05 


-.0021(.0023) - 


0227(.0021) 


0929(.0139) 


.0053( 


.0003) 


.0050( 


.0009) 


.2073( 


.0182) 


3.12 


0.03 


23 


.50 


95 


.05 


.0015(.0024) - 


0237(.0026) 


1005(.0136) 


.0060( 


.0004) 


.0075( 


.0024) 


.2022( 


.0198) 


3.12 


0.03 


34 


.50 


95 


.05 


.0020(.0025) - 


0222(.0022) 


1045(.0137) 


.0063( 


.0006) 


.0054( 


.0013) 


.2043( 


.0170) 


3.12 


0.03 


13|2 


.50 


95 


.05 


-.0084(.0024) - 


0190(.0017) 


0459(.0127) 


.0062( 


.0004) 


.0032( 


.0005) 


.1695( 


.0131) 


3.12 


0.03 


24 3 


.50 


95 


.05 


-.0131(.0023) - 


0205(.0021) 


0332(.0137) 


.0056( 


.0003) 


.0048( 


.0011) 


.1940( 


.0177) 


3.12 


0.03 


14|23 


.50 


95 


.05 


-.1324(.0022) - 


0245(.0021) - 


0670(.0137) 


.0227( 


.0008) 


.0053( 


.0013) 


.1980( 


.0145) 


3.12 


0.03 



Table 1: Estimated relative bias and relative MSE with estimated standard errors for four dimensional D-vine SCAR copulas 
with common bivariatc copula family and common time-varying structure for all pair copulas. 



the copula on the second tree. Again we report estimated relative bias and relative MSE together with 
the signal-to-noisc and stationary variance of the latent AR(1) process in Table [31 As before, the results 
are based on 1000 simulated data sets of length 1000. 

The results in Table [3] show only a moderate influence of different time- varying structures or different 
pair copula families. It is notable that lower persistence leads to more precise estimation of /i and (\) and 
less precision for a . but that effect can mainly be explained by the changed signal-to-noise ratio. Again 
parameters of the higher tree are less well estimated compared to the ones in the first tree. 

Conclusions from the simulation study 

Overall we saw that the estimation becomes worse on higher trees and that this effect is stronger 
when we allow for time variation. The relative imprecision also increases with lower overall dependence. 
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Ind 


True values 




Average relative Bias 






Average relative MSE 




sn 


avar 








4' 


(T 


/J 


<P 




a 


M 








a 






Copula families: Gaussian Gaussian Gaussian Gaussian Gaussian Gaussian 
















12 


.50 


.95 


.15 


.0395( 


.0083) -.0109(.0009) 


-.0827 


.0035) 


.0689{ 


.0075) 


.0009( 


.0001) 


.0192 


.0011) 


1.04 


0.23 


23 


.50 


.95 


.15 


.0291( 


.0070) -.0127(.0009) 


-.0721 


.0033) 


.0493{ 


.0034) 


.0010( 


.0001) 


.0158 


.0007) 


1.04 


0.23 


34 


.50 


.95 


.15 


.0263( 


.0078) -.0113(.0008) 


-.0763 


.0034) 


.0598{ 


.0055) 


.0008( 


.0001) 


.0171 


.0009) 


1.04 


0.23 


13|2 


.50 


.00 


.00 


-.0156( 


.0012) -(- ) 




- ) 


.0016{ 


.0001) 








-) 


Inf 





24 3 


.30 


.00 


.00 


.0416( 


.0022) -(- ) 




- ) 


.0067{ 


.0003) 




_\ 




-) 


Inf 





14|23 


.20 


.00 


.00 


.2185( 


.0039) -(- ) 




- ) 


.0624{ 


.0017) 








-) 


Inf 





Copula families : Clayton 


Clayton 


Clayton Gaussian Gaussian 


Gaussian 
















12 


.50 


.95 


.15 


.0781( 


.0088) -.0077(.0011) 


-.0403 


.0041) 


.0760{ 


.0095) 


.0011( 


.0002) 


.0166 


.0011) 


1.04 


0.23 


23 


.50 


.95 


.15 


.0712( 


.0088) -.0110(.0012) 


-.0326 


.0041) 


.0751( 


.0115) 


.0014( 


.0002) 


.0162 


.0012) 


1.04 


0.23 


34 


.50 


.95 


.15 


.0544( 


.0091) -.0091(.0011) 


-.0347 


.0045) 


.0777{ 


.0103) 


.0012( 


.0002) 


.0192 


.0014) 


1.04 


0.23 


13|2 


.50 


.00 


.00 


-.0152( 


.0012) -(- ) 




- ) 


.0015{ 


.0001) 




_) 




-) 


Inf 





24 3 


.30 


.00 


.00 


.0424( 


.0023) -(- ) 




- ) 


.0067( 


.0003) 




_) 




-) 


Inf 





14|23 


.20 


.00 


.00 


.2237( 


.0043) -(- ) 




- ) 


.0671( 


.0020) 




\ 




-) 


Inf 





Copula families : Gumbel Gumbel 


Gumbel Gaussian Gaussian 


Gaussian 
















12 


.50 


.95 


.15 


.0282( 


.0063) -.0043(.0007) 


-.0259 


.0038) 


.0404{ 


.0022) 


.0005( 


.0001) 


.0153 


.0007) 


1.04 


0.23 


23 


.50 


.95 


.15 


.0243( 


.0065) -.0055(.0007) 


-.0193 


.0040) 


.0421{ 


.0023) 


.0005( 


.0001) 


.0161 


.0010) 


1.04 


0.23 


34 


.50 


.95 


.15 


.0114( 


.0069) -.0039(.0006) 


-.0251 


.0039) 


.0473{ 


.0055) 


.0004( 


.0000) 


.0157 


.0013) 


1.04 


0.23 


13|2 


.50 


.00 


.00 


-.0070( 


.0011) -(- ) 




- ) 


.0013{ 


.0001) 








-) 


Inf 





24 3 


.30 


.00 


.00 


.0500( 


.0022) -(- ) 




- ) 


.0074{ 


.0003) 




_\ 




-) 


Inf 





14|23 


.20 


.00 


.00 


.2346( 


.0038) -(- ) 




- ) 


.0696{ 


.0018) 




_\ 




-) 


Inf 





Copula families: Gaussian Gaussian Gaussian Gaussian 


Gaussian Gaussian 
















12 


.50 


.95 


.05 


-.0012( 


.0025) -.0168(.0020) 


.0465 


.0112) 


.0064{ 


.0007) 


.0045( 


.0024) 


.1331 


.0105) 


3.12 


0.03 


23 


.50 


.95 


.05 


.0032( 


.0023) -.0179(.0017) 


.0575 


.0108) 


.0054{ 


.0003) 


.0034( 


.0009) 


.1266 


.0087) 


3.12 


0.03 


34 


.50 


.95 


.05 


.0037( 


.0022) -.0160(.0018) 


.0455 


.0106) 


.0052{ 


.0002) 


.0036( 


.0014) 


.1191 


.0079) 


3.12 


0.03 


13|2 


.50 


.00 


.00 


.0059( 


.0012) -(- ) 




- ) 


.0015{ 


.0001) 


_l 


_) 


_ 


-) 


Inf 





24|3 


.30 


.00 


.00 


.0295( 


.0021) -(- ) 




- ) 


.0053{ 


.0002) 




_) 




-) 


Inf 





14|23 


.20 


.00 


.00 


.2502( 


.0033) -(- ) 




- ) 


.0740{ 


.0017) 




_\ 




-) 


Inf 





Copula famil 


es: Clayton 


Clayton 


Clayton Gaussian Gaussian Gaussian 
















12 


.50 


.95 


.05 


-.0014( 


.0024) -.0191(.0023) 


.0631 


.0107) 


.0060{ 


.0008) 


.0060( 


.0023) 


.1233 


.0132) 


3.12 


0.03 


23 


.50 


.95 


.05 


.0020( 


.0022) -.0155(.0013) 


.0573 


.0099) 


.0051{ 


.0002) 


.0019( 


.0002) 


.1061 


.0068) 


3.12 


0.03 


34 


.50 


.95 


.05 


.0041( 


.0022) -.0135(.0012) 


.0461 


.0096) 


.0052{ 


.0002) 


.0016( 


.0002) 


.0994 


.0070) 


3.12 


0.03 


13 


2 


.50 


.00 


.00 


.0061( 


.0012) -(- ) 




- ) 


.0015{ 


.0001) 








-) 


Inf 





24 


3 


.30 


.00 


.00 


.0284( 


.0021) -(- ) 




- ) 


.0052{ 


.0002) 








-) 


Inf 





14123 


.20 


.00 


.00 


.2482( 


.0034) -(- ) 




- ) 


.0740{ 


.0018) 








-) 


Inf 





Copula 


families : Gumbel Gumbel 


Gumbel Gaussian Gaussian 


Gaussian 
















12 


.50 


.95 


.05 


-.0010( 


.0026) -.0211(.0022) 


.0902 


.0136) 


.0069{ 


.0007) 


.0056( 


.0023) 


.2006 


.0163) 


3.12 


0.03 


23 


.50 


.95 


.05 


.0036( 


.0023) -.0247(.0026) 


.1116 


.0135) 


.0054{ 


.0002) 


.0079( 


.0031) 


.2035 


.0144) 


3.12 


0.03 


34 


.50 


.95 


.05 


.0044( 


.0024) -.0223(.0023) 


.0964 


.0133) 


.0058{ 


.0004) 


.0061( 


.0020) 


.1935 


.0161) 


3.12 


0.03 


13 


2 


.50 


.00 


.00 


.0042( 


.0012) -(- ) 




- ) 


.0014{ 


.0001) 








-) 


Inf 





24 


3 


.30 


.00 


.00 


.0262( 


.0021) -(- ) 




- ) 


.0051{ 


.0002) 








-) 


Inf 





14|23 


.20 


.00 


.00 


.2472( 


.0034) -(- ) 




- ) 


.0730{ 


.0018) 








-) 


Inf 






Table 2: Estimated relative bias and relative MSE with estimated standard errors for four dimensional D-vine SCAR copulas 
with common bivariatc copula family and common time- varying structure for all first tree pair copulas. 



index True values 

(T 


relative bias 
/i cr 


relative MSE 
/I cr 


sn avar 


Copula families: Gaussian 
12 .50 .85 .15 
23 .50 .95 .15 
13|2 .50 .85 .15 


Gaussian Gaussian 

-.0064(.0023) .0035(.0015) -.0864(.0045) 
.0403(.0073) -.0117(.0009) -.0774(.0035) 
-.0379(.0021) -.0022(.0020) -.1801(.0049) 


.0057(.0003) .0023(.0001) .0280(.0013) 
.0567(.0050) .0009(.0001) .0190(.0010) 
.0059(.0003) .0041(.0003) .0577(.0020) 


1.76 0.08 
1.04 0.23 
1.76 0.08 


Copula families : Clayton 
12 .50 .85 .15 
23 .50 .95 .15 
13|2 .50 .85 .15 


Clayton Clayton 
-.0026(.0023) -.0014(.0014) -.0541(.0044) 
.0696(.0075) -.0091(.0010) -.0335(.0041) 
-.1193(.0018) -.0067(.0022) -.2900(.0049) 


.0053(.0002) .0021(.0001) .0221(.0010) 
.0608(.0079) .0011(.0002) .0179(.0015) 
.0176(.0005) .0048(.0003) .1085(.0029) 


1.76 0.08 
1.04 0.23 
1.76 0.08 


Copula families: Gumbel 
12 .50 .85 .15 
23 .50 .95 .15 
13|2 .50 .85 .15 


Gumbel Gumbel 
-.0034(.0024) -.0053(.0017) -.0410(.0054) 
.0245(.0061) -.0038(.0006) -.0254(.0037) 
-.0588(.0020) .0036(.0020) -.2254(.0054) 


.0058(.0003) .0030(.0002) .0312(.0014) 
.0393(.0021) .0004(.0001) .0145(.0008) 
.0078(.0003) .0043(.0003) .0806(.0026) 


1.76 0.08 
1.04 0.23 
1.76 0.08 


Copula families : Clayton 
12 .50 .85 .15 
23 .50 .95 .15 
13|2 .50 .85 .15 


Gumbel Gaussian 
-.0018(.0023) -.0013(.0014) -.0545(.0043) 
.0290(.0064) -.0031(.0006) -.0270(.0037) 
-.0372(.0020) -.0038(.0022) -.1854(.0050) 


.0053(.0002) .0021(.0001) .0222(.0010) 
.0430(.0025) .0004(.0001) .0146(.0007) 
.0057(.0002) .0049(.0005) .0606(.0020) 


1.76 0.08 
1.04 0.23 
1.76 0.08 



Table 3: Estimated relative bias and relative MSE together with estimated standard errors in parentheses for selected 
senarios of three dimensional SCAR D-vine copulas assuming different time-varying structures and common or different 
pair copula families. 
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For our application we can conclude that it may be sensible to restrict the time variation to the first few 
trees and restrict the copula parameter to be constant beyond. Furthermore, placing the pairs with the 
largest dependence on the first tree as is commonly done for D-vine copula models is expected to provide 
the most precise estimates for the time variation in dependence, which is what we arc mainly interested 
in here. 

5. Application 

In this section we provide an empirical illustration of the D-vine SCAR model. The dataset we consider 
are daily returns from 29 stocks listed in the DAX30 index during the period from the 1"* of January 2007 
to the 14*'* of June 2011. giving a total of 1124 observations. A list of the included companies is given 
in the Appendix. We decided for this dataset to find a balance between demonstrating the possibility of 
high dimensional modeling and the ability to still present the main results. Nevertheless, one could in 
principle consider much larger dimensions, which we leave for future research. 

The conditional mean of the returns was modeled using an ARMA(1,1) model for all stock separately. 
The Ljung-Box statistics for the residuals revealed no significant remaining autocorrelation. For the 
conditional variance we considered GARCH(1,1) models with Student t innovations. Although one may 
argue that GARCH models that allow for the leverage effect such as the GJR-GARCH arc appropriate 
for many individuals stocks, preliminary results suggested that the results for the dependence modeling 
arc not affected by this choice. Results for the marginal models are not reported here for brevity but are 
available upon request. 

Next, we estimated the D-vine SCAR model on the transformed standardized residuals. The ordering 
of the variables was done by maximizing the overall pairwisc dependence measured by Kendall's tan. 
In particular, first choose the pair of variables with the highest empirical Kendall's t. Second connect 
the next variable which has highest pairwisc Kendall's r with one of the previously chosen variables 
and proceed in a similar fashion until all variables are connected. This is the common strategy for D- 
vinc copulas and is also motivated by our simulation results. In particular, we expect to capture the 
overall time variation of the dependence as good as possible with this choice, as it turns out that time 
variation is most relevant on the first tree. For each bivariate (conditional) copula model we then face 
two important choices, namely whether dependence is time- varying or static, and which copula family to 
use. We automatically select the model by first estimating time-varying and constant copulas from the 
following families: Gumbel (G), survival Gumbel (SG), Clayton (C), survival Clayton (SC), Normal (N) 
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and the independence copula (ijj. We then select the best fitting copula family from these 11 candidate 
models by the Bayesian information criterion (BIG). Given the size and complexity of our model, as 
well as the difficulty to estimate parameters precisely on higher trees, we decided to rely on the BIG 
to find more parsimonious model specifications and to minimize the estimation error. Additionally, we 
considered the restriction of allowing potential time variation only on a limited (small) number of trees. 
We considered this restriction from 1 to 12 trees, but we report only a subset of these models since results 
are identical or at least very similar for many of those cases. Specifically, it turned out that making this 
restriction beyond the 6*'* tree is irrelevant, since there is no evidence of time variation on higher trees. 
A further possible restriction that may be made is to truncate the vine beyond a certain tree, meaning 
that all conditional copulas are set to the independence copula. In the current application we did not 
make this restriction because the independence copula is included in the set of admissible models and 
the automatic selection by the BIG in practice leads to a truncation. Nevertheless, for large dimensional 
applications truncation of the vine should definitely be considered. 





MO 


Ml 


M2 


M3 


M4 


M6 


time- varying 




1 


1-2 


1-3 


1-4 


1-6 


time-constant 


1-27 


2-27 


3-27 


4-27 


5-27 


7-27 


up to tree 1 


-11491 


-12177 


-12177 


-12177 


-12177 


-12177 


up to tree 2 


-13657 


-14313 


-14346 


-14346 


-14346 


-14346 


up to tree 3 


-14893 


-15482 


-15517 


-15525 


-15525 


-15525 


up to tree 4 


-15653 


-16207 


-16239 


-16265 


-16265 


-16265 


up to tree 5 


-16068 


-16538 


-16563 


-16583 


-16584 


-16584 


up to tree 6 


-16326 


-16774 


-16787 


-16796 


-16802 


-16804 


up to tree 7 


-16466 


-16882 


-16892 


-16892 


-16909 


-16914 


up to tree 20 


-17358 


-17621 


-17599 


-17580 


-17585 


-17594 


total 


-17366 


-17626 


-17604 


-17587 


-17593 


-17603 



Table 4: Partial BIG (up to trees 1-7 and 20) and total BIG values for Models with no (MO), first tree (Ml), first and 
second (M2) , first three (M3), first four (M4) and first six trees (M6) time- varying 



The results of our estimation presented in Tables H] and [5] report the estimation results when allowing 
time variation in no, in trees 1 to 4 and in tree 6, respectively. For trees 5 and 7-28 no time-varying 
dependence was found. Table [S] reports the number of time- varying copulas found on each tree and how 
often each copula family was selected. In addition, the number of parameters estimated in each tree 
is given. In Table |3] partial and total BIG values are reported. We note that time variation is very 
important, when modeling the dependence in the first tree. Here 27 out of 28 copulas are chosen to be 
time-varying. Furthermore, the vast majority of the time-varying copulas is selected to be the Gaussian 



''Obviously, for the independence copula no parameter needs to be estimated. 
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copula, which is in line with the fi ndings 
given in 



Hafner and Manned (|2011D . An explanation for this finding is 



Manner and Segera (1201 II ). who show that the Gaussian copula with random correlations has 
much larger dependence in the tails than the static Gaussian copula. This is in line with the stylized 
fact that financial returns arc characterized by tail dependence. Beyond the first tree, however, only very 
few pairs have time-varying dependence justifying the restrictions. Among the copulas with constant 
parameters all families arc chosen, but Gaussian, survival Gumbel and independence copulas arc selected 
more often than the other families. Especially on higher trees the independence copula dominates, 
indicating that a truncation after a certain level (say 15-20) would come almost without any costs in 
terms of model fit. The overall fit of the models measured by the BIG for the total model turns out to be 
best when allowing time variation only for on the first level. This can be explained by the strong penalty 
on additional parameters by the BIG and the better fit on the higher trees. 

Finally, we are interested in estimating the path of the pairwise dependence parameters. Smoothed 
estimates for the time-varying Kendall's r based on the model with time-variation only on the first 
tree are presented in Figures [5] and |3] for the pairwise dependence of the five companies that have the 
largest weight in the DAX index. The companies are (with their corresponding node in the D-vine in 
parenthesis) AUianz (13), Bayer (9), E.ON (12), SAP (27) and Siemens (19). This choice covers the 
interesting situations that we have companies that are neighboring in the vine, i.e. their dependence 
is modeled directly, that companies are close in the vine and that they lie quite far apart from each 
other. In the latter situation the implied pairwise time- varying dependence is computed conditional on 
the dependence between all variables in between. Note that in this case the dependence parameter at 
each point in time has to be computed using Monte Garlo simulation, for which we used 400 Monte Garlo 
replications. The dynamics of the dependence parameter arc clearly visible for all situations and we 
have strong evidence of dependence changing over time. For example, the dependence parameters have 
decreased strongly in 2009 for most pairs. It also stands out that the dependence parameters involving 
the company Bayer have much less pronounced movements in dependence, which is mostly negative or 
close to independence. 



6. Conclusions and further research 

In light of the current financial crisis including the discussion of understanding systemic risk the need 
to understand time-varying effects not only within individual financial products but also among groups 
of financial variables has been increasing. The developed D-Vine SCAR models are aiming to fill this 
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Time Varying Kendall's tau between Alllanz and Bayer 



Time Varying Kendall's tau between Alllanz and E-ON 




— I 1 1 \ 1 — I 1 \ 1 1 — 

2007 2008 2009 2010 2011 2007 2008 2009 2010 2011 

Time Varying Kendall's tau between Bayer and E-ON Time Varying Kendall's tau between Bayer and SAP 



I 1 \ 1 1 

2007 2008 2009 2010 2011 2007 2008 2009 2010 2011 

Figure 2: Smoothed estimates for time varying Kendall's r 

demand. From a statistician's point of view such a model is demanding. First a very flexible multivariate 
dependency model is required such as the class of vine copulas and secondly an appropriate model for 
the time dependency of the copula parameters has to found. For parameter interpretability we chose to 
follow a parameter driven instead of an observation driven approach and thus allowed for a stochastic 
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Time Varying Kendaii's tau between Bayer and Siemens 



Time Varying Kendaii's tau between E-ON and SAP 




Time Varying Kendaii's tau between E-ON and Siemens 



Time Varying Kendaii's tau between SAP and Siemens 





Figure 3: Smoothed estimates for time varying Kendall's r 



autoregressive structure for the copula parameters to model time-varying dependence among a large 
group of variables. 

While this approach leads to a relatively straightforward model formulation the development of effi- 
cient estimation procedures is much more difficult. Especially in high dimensions maximum likelihood is 
infeasible, since it would require the maximization over integrals of size equal to the data length. In the 
application the data length was 1125. These integrals occur since we need to integrate over the latent 
variable process to express the joint likelihood. This problem already occurs when we consider bivariat e 



SCAR models. One solution to this is to use efficient importanc e samp ling (jRichard and Zhang 



Aasetal 



20071) 



(|2009l ) for multivariate copulas allows 



In addition, the pair copula construction approach of 
to express the likelihood in bivariate copula terms in addition to a sequential formulation over the vine 
tree structure. This makes it feasible to develop and implement efficient importance sampling for the 



D-vine SCAR model. In a simulation study we validated our estimation procedure and the application 
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to joint dependency modeling of 29 stocks in the DAX showed that time- varying dependence structures 
can be observed. One interesting feature of this apphcation is that nonnormal pair copulas with constant 
parameters were replaced by normal pair copulas when time- varying copula parameters are allowed. This 
was observed when the strength of the dependence was moderate to large. 

In this paper we followed some first approaches to model selection. We first restricted to a known 
dependency structure given by a D-vine, but allowed the copula family of each pair copula to be chosen 
among a prespecified class of copula families in addition to the choice if a pair copula has time- varying 
parameters or not. For this we used BIG, however more sophisticated criteria might be necessary. We 
also restricted the use of time- varying pair copul a parameters to a prespe cified number of top trees. Here 



the approach of truncated vines as developed in 



Brechmann et al. 



(|201l[ ) might be a good starting point 



to choose this number in a data driven manner. As already mentioned it is feasible to extend the class 
of D-vine SCAR models to include R-vines as copula models. More research is also needed to investigate 
the effects of time- varying copula models on economic entities such as portfolio returns and value at risk. 
Finally the model uncertainty introduced by assuming the marginal par ameter estimates as true ones in 
the two-step approach has to assessed. However the simulation results of |] 
valid for a copula model with time-varying parameters. 



Kim et al 



(|2007l ) might remain 
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MO: all trees time-constant 



Tree 


time- varying 




time-constant 




# par 




N 


SG 


I 


N 


C 


G 


SC 


SG 




1 









7 





4 





17 


28 


2 









10 





12 





5 


27 


3 






2 


14 





7 





3 


24 


4 






2 


15 





4 


1 


3 


23 


5 






4 


15 





4 





1 


20 


Ml: tree 1 time- varying, trees 2-27 time-constant 




1 


23 


4 

















1 


82 


2 









4 





8 





15 


27 


3 






1 


9 





5 


1 


10 


25 


4 






2 


13 





7 





3 


23 


5 






5 


16 





2 


1 





19 


M2: trees 1-2 time-varying 


, trees 3-27 time-constant 


1 


23 


4 

















1 


82 


2 


1 


2 





4 





8 





12 


33 


3 






1 


8 





5 


1 


11 


25 


4 






2 


14 





6 





3 


23 


5 






5 


15 





2 


1 


1 


19 


M3: trees 1-3 time-varying, trees 4-27 time-constant 


1 


23 


4 

















1 


82 


2 


1 


2 





4 





8 





12 


33 


3 


3 





1 


7 





5 


1 


9 


31 


4 






2 


10 





7 





6 


23 


5 






5 


14 


1 


2 


1 


1 


19 


M4: trees 1-4 time-varying, trees 5-27 time-constant 


1 


23 


4 

















1 


82 


2 


1 


2 





4 





8 





12 


33 


3 


3 





1 


7 





5 


1 


9 


31 


4 


1 





2 


9 





7 





6 


25 


5 






5 


14 


1 


2 


1 


1 


19 


M6: trees 1-6 time-varying 


, trees 7-27 time-constant 


1 


23 


4 

















1 


82 


2 


1 


2 





4 





8 





12 


33 


3 


3 





1 


7 





5 


1 


9 


31 


4 


1 





2 


9 





7 





6 


25 


5 








5 


14 


1 


2 


1 


1 


19 


6 


1 





11 


6 


1 


1 





3 


14 


7 






8 


4 


1 


3 





6 


14 



Table 5: Number of bivariate copula families used with time- varying and time-constant parameters and total number of 
parameters in trees 1-5 for models allowing for increasing number of trees with time-varying parameters 
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Appendix A. Stocks in the DAX and their ordering 



Company 


Ticker Symbol 


node in D-vine 


Adidas 


ADS 


25 


AUianz 


ALV 


13 


BASF 


BAS 


- 


Bayer 


BAYN 


9 


Beiersdorf 


BEX 


6 


BMW 


BMW 


18 


Commerzbank 


CBK 


15 


Daimler 


DAI 


17 


Deutsche Bank 


DBK 


14 


Deutsche Brse 


DBl 


24 


Lufthansa 


LHA 


16 


Deutsche Post 


DPW 


22 


Deutsche Telekom 


DTE 


10 


E.ON 


EOAN 


12 


Fresenius 


PRE 


3 


Fresenius Medical Care 


FME 


2 


HeidelbergCement 


HEI 


26 


Henkel 


HENS 


7 


Infineon Technologies 


IPX 


28 


K+S 


SDF 


4 


Linde 


LIN 


8 


MAN 


MAN 


21 


Merck 


MRK 


1 


Metro 


MEO 


5 


Munich Re 


MUV2 


23 


RWE 


RWE 


11 


SAP 


SAP 


27 


Siemens 


SIE 


19 


ThyssenKrupp 


TKA 


20 


Volkswagen Group 


VOWS 


29 



Table A. 6: DAX companies and their node position in the selected D-vine for the sample period from the f^' of January 
2007 until the 14*'' of June 2011. The company BASF(BAS) was dropped from the analysis, because data was not available 
for the complete sample period. 
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